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ABSTRACT 

Context. X-ray spectra of AGN often contain signatures indicative of absorption in multiple layers of gas whose ionization-state and covering 
fraction may vary with time. It has been unclear to date how much of the observed X-ray spectral and timing behavior in AGN can be attributed 
to variations in absorption, versus variations in the strengths of emission or reflection components. Diagnostics of the inner regions of AGN 
cannot be reliably performed until the origin of observed effects is understood. 
Aims. We investigate the role of the X-ray absorbers in the Seyfert 1 galaxy NGC 3516 

Methods. Time-averaged and flux-selected spectroscopy is used to examine the behavior of NGC 3516 observed in Chandra HETG and XMM 
data from Oct 2006. 

Results. New H-like and He-like emission and absorption features discovered in the Fe K regime reveal a previously unknown zone of 
circumnuclear gas in NGC 3516 with log £ ~ 4.3 and column density ~ 10 23 cnr 2 . A lower-ionization layer with log £ ~ 2 and of similar 
column density is confirmed from previous observations, this layer has a covering fraction around 50%, and changes in covering provide a 
simple explanation of a deep dip in the light curve that we interpret as an eclipse of the continuum due to passage of a cloud across the sight 
line within half a day. These inner zones of absorbing gas are detected to have outflow velocities in the range 1000 - 2000 km s~', this, and 
constraints on radial location are consistent with an origin as part of a disk wind in NGC 3516. 
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1. Introduction 



The process of accretion onto a supermassive black hole feeds 
the so-called 'active' subset of galaxies. The centers of galax- 
ies are regions where material must exist in the strong gravity 
regime and in active galaxies there is also abundant radiation 
emerging in the X-ray bandpass from that region. The inter- 
action of nuclear X-rays with circumnuclear gas allows us a 
direct probe of the extreme conditions close to an event hori- 
zon. It is widely thought that the large amount of angular mo- 
mentum that accreting material must have means this material 
forms an 'accretion disk' around the black hole, extending to 
within a few gravitational radii (the innermost stable orbit de- 
pending on the angular momentum of the black hole). Other 
significant components likely exist in the nuclear environs in- 
cluding a corona of high energy particles sandwiching the disk, 
a wind of material ablated from the disk surface and clouds of 



Intrinsically narrow emission lines arising from the surface 
of the accretion disk or in any of the gas zones very close to 
the black hole will be distorted by a combination of Doppler 
and relativistic effects w hose details can yi eld conditions and 
kinematics of the gas (see lFabian etai]2000l for a review of ef- 
fects). X-ray spectroscopy provides a good probe of gas very 
close to the black hole as Fe Ka is a strong emission line, emit- 
ted in the X-ray band via fluorescence or recombination pro- 
cesses between 6.4 - 7 keV (depending on the ionization-state 
of the gas); this is a strong line owing to the high abundance 
and fluorescence yield of Fe. The high ionization potentials of 
the highest ionization states of iron mean that this line can be 
produced in regions where the gas is extremely ionized. 

Fe Ka emission is commonly observed in AGN and much 
has been made of the diagnostic potential of the line, with many 
authors deriving accretion di sk parameters from the line profile 



gas a t various radial locations. ICrenshaw. Kraemer & George 
(2003) review the details of X-ray reprocessors along with esti- 
mates of the locations and conditions in various distinct zones. 



with important implications (Tan aka et all 1995; Na ndra et all 



1997b . However, spectroscopy at the modest energy resolution 
currently available in the X-ray band has been plagued with 
ambiguity; complex absorption from ionized gas can produce 
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curvature in the 2-6 keV band that looks similar to that of an 
Fe emission line produced close to the event horizon of the ac- 
cretion disk. Recent spectra from XMM have shown absorp- 
tion lines from highly ionized ions in several AGN, allow- 
ing construction of detailed absorption models in some cases 
and revealing hith erto unknown zones of absorbing gas (e.g. 
Reeves et al.l2004 . In th e broad-line S eyfert 1 galaxy (BLSyl) 
NGC 3516 (z=0.008836: lKeele II 1996b . absorption lines from a 
broad range of ions showed three distinct zones of gas covering 
the active nucleus in that case and producing muc h of the spec 



X-ray observations of NGC 35 16 also provided the first ev- 
idence for 'transient' narrow emission lin es, most eas i ly ex- 



plained as Doppler shifted Fe em ission (iTurner. et all 12002 



tral c urvature and variability in the X-ray band ([Turner et al. 



2005b . The classic U versi on of ionization parameter was uti- 
lized bv lTurner et al. (2005) Here we use the £ form of ioniza- 
tion parameter where £ = L/nr 2 . L is the 1-1000 Rydberg lu- 
minosity, n the gas density and r the absorber-source distance. 
We quote £ in units of ergcms -1 throughout and xstar tables 
are based on runs with micro-turbulence <x = 200 km s _1 unless 
otherwise noted. Ionization states found for the the gas zones 
present during 2001 observations can be converted into £ val- 
ues assuming an illuminating intrinsic X-ray continuum r ~ 2 
with a -x ~ 1.5. Thus the 2001 results correspond to: 



- a cool 'UV absorber' with log %\ 
10 21 cnT 2 

- an outflowing zone with log £2 ~ 3, Nh 
outflow velocity ~ 1 100 km s -1 

- a 'heavy' absorber with log £3 ~ 2.5, Nh ~ 10 23 crrT 2 
ering ~ 50% of the continuum. 



-0.5, N H ~ 6 x 
10 22 cirT 2 and 
cov- 



Turn er et all (|2005) found the spectral variability observed 



within 2001 to be consistent with the absorbing gas layers re- 
sp onding in ionization-state as the continuum flux varied. 



Mar kowitz et al. d2007b presented a 2005 Suzaku observa- 



tion where two zones of gas were detected: the primary ab- 
sorber is a layer of low-ionization log £ ~ 0.3 and Nh ~ 
6 x 10 22 cirT 2 covering ~ 96 - 100% of the source, while ev- 
idence is also found for a zone having log £ ~ 3 and Nh ~ 
4 x 10 22 cirT 2 . It is unclear which of these zones correspond to 
particular layers isolated in 2001 data. Suzaku data were con- 
sistent with the presence of the UV-absorber detected during 
2001, but if we want to associate the primary absorber detected 
by Suzaku with the outflowing layer detected in 2001, then the 
gas has undergone a major change in ionization-state between 
the two epochs; if we compare the Suzaku primary absorber 
with the 'heavy' layer detected during 2001, then that layer 
would have had to have changed Nh and covering fraction. 
Whatever the correspondence between modeled zones Suzaku 
data clearly show NGC 3516 to have much larger absorption 
by re latively low-ionization gas compared to any of the 2001 
data dMarkowitz et alll2007l) . During 2006 Suzaku found NGC 
3516 to have hard-band flux (above ~ 6 keV) comparable to 
that observed during April 2001 and thus the change in absorp- 
tion seems likely to be due to the passage of a relatively cool 
cloud into the sight-line, as s een in X-ray observa tions of some 
other AGN (e.g. N GC 4151 iKraemer et al1l2005l NGC 1365 
Risal iti et alll2007l) . At this point one could only conclude that 
the inner zones of the absorber in NGC 3516 change dramati- 
cally over time in all senses. 



Iwasawa. Miniutti & Fabian , 2004b . With all of these new dis- 
coveries our picture of the contributions from various reproces- 
sors is rapidly evolving. 

In this paper we present data from a coordinated observa- 
tion of NGC 35 16 utilizing Chandra and XMM data, exploring 
further the innermost absorbing gas in this source. 

2. The Observations 

An XMM-Newton (hereafter XMM) observation of NGC 3516 
was made covering 2006 Oct 6-13 (OBSIDs 04012010401, 
501, 601, 1001). Chandra observations (OBSIDs 7281,7282, 
8450, 8451, 8452) were interwoven with the XMM exposures, 
to provide complete coverage of the source during that period. 
Fig.Q] shows the light curve from the two satellites, illustrat- 
ing the coverage obtained and the source behavior across the 
observation. 

EPIC data utilized the thinl filter and pn 
PrimeSmallWindow mode. Data were processed using 
SAS v7.0 using standard criteria and removing periods of high 
background. Instrument patterns 0-4 were selected for the pn 
and here our analysis describes only the pn data, as these offer 
the best S/N ratio of all EPIC data, and were completely free 
of photon pileup effects. Data were extracted from a circular 
cell of radius 40" centered on the source, and background data 
were taken from a source-free region of the same pn chip. 
Periods where the full-band pn count rate exceeded 2 count s _1 
in the background region were excluded from the analysis. The 
combined screening yielded a total effective EPIC pn exposure 
of 155 ks over an observational baseline of about one week. 

NGC 3516 gave a mean pn count rate ~ 4.91 +0.007 
counts 1 in the 2-10 keV band. The mean background level 
in the screened data was ~ 1 % of the mean source rate in this 
band. As evident in Fig.[T] the source ranges across a 2-10 keV 
flux F ~ 1.4 - 5.4 x 10~ n erg crrr 2 s~' during the course of the 
observations. First-order RGS data yielded 0.450+0.001 and 
0.389+0.001 counts 1 in RGS 1 and 2, respectively. The RGS 
background level was ~ 12% of the total count rate. 

The Chandra HETG exposure yielded a 2-7 keV count rate 
~ 0.021 ±0.0005 count s _I in the summed HEG first order spec- 
trum. Data were reduced using CIAOv3.4 and CALDB v3.4.0 
and following standard procedures for extraction of HETG 
spectra, with the exception that we utilized a narrower extrac- 
tion strip than the tgextract default. The default processing 
criterion cuts off the HEG data above 8 keV, yet some of the 
brighter and/or well-exposed sources still have good HEG data 
above 8 keV, including this well exposed observation of NGC 
3516. This reason for the default processing cut-off is that the 
overlap of the MEG and HEG strips depends on the extraction 
strip widths, and if the latter are too large, a larger intersection 
of the MEG and HEG strips results, cutting off the HEG data 
prematurely. Specifically, we used width_factor_hetg=2Q 
in the tool tg_create_mask, instead of the default value of 
35. As the HETG spectra have a very low signal-to-noise ra- 
tio it was necessary to coadd the positive and negative first- 
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Fig. 1. XMM pn data in 500 s bins, and summed HEG and 
MEG data from all orders over 1-6 keV in 1500 s bins. The 
light curve shows the coverage of NGC 3516 with each satel- 
lite. The XMM data were time-split for spectral analysis, and 
the red points denote the time selection defining the high-flux 
state; the black points show the medium or mid-flux state, the 
green points show the dip event, comprising the low-flux state 

order spectra, and coadd all five OBSIDs, to create a high 
quality summed first-order HEG and MEG spectra for fitting. 
After such co-addition, we binned the spectra to 4096 channels, 
which yielded > 20 counts per bin in the fitted range. 

3. Results 

The light curve from the combined observation shows the 
source to exhibit significant flux variability (Fig.[TJi, with a no- 
table flare early in the observation, and a drop into a 'dip' state 
later in the observation (Fig.|2j». Comparison of spectra from 
the different states show a flattening as the source flux drops 
(Fig. O as often observed in sources of this class. Comparison 
of the new spectra with 2001 XMM pn and 2005 Suzaku XIS 
data illustrates the relatively high flux level at which the source 
was observed during 2006 (Fig.|4]i. 

Inspection of the detailed profile of the light curve during 
the dip reveals interesting structure (Fig.[2]i, with the source 
flux levelling off for a period lasting ~ 10 ks. Intriguingly, the 
initial shallow dip egress has a duration of 20 ks, but breaks to 
a sharp rise occuring over just 6 ks. 

We first analyzed the time-averaged pn spectrum. Fig. [5] 
shows a fit to the pn data in the Fe K-band with 2006 data 
shown in red. The most obvious result is the presence of two 
deep absorption lines in the 6.5-7 kev regime, naturally iden- 
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Fig. 2. A close-up view of the dip event from Fig. 1 . 



tified with He-like and H-like species of Fe. Fitting the line 
pair with simple gaussian profiles we find an equivalent width 



100+JI eV for the line at E 



6.989 keV, and 45+'° eV for 



the line at E — 6.697 keV, both measured against the total 'lo- 
cal continuum' . To obtain the most accurate line energies and 
widths we turned to the contemporaneous HETG data. 

Fig. |6] shows the HEG spectrum of NGC 3516 from 2006. 
The deep absorption lines from Fe are confirmed, as are narrow 
emission lines in the same band. Of great interest is the detec- 
tion of an emission line redward of the H-like Fe absorption 
feature. 

Fitting the emission and absorption line 'pair' close to 7 
keV yields energies close to those of rest-frame H-like Fe tran- 
sitions; the red-shifted emission and blue-shifted absorption of 
these components constitutes a profile somewhat akin to a clas- 
sic 'P Cygni' profile, indicative of an origin in an outflowing 
layer of gas, although in this case the emission line is offset 
from the expected rest energy, as well as the absorption line. 
A redshifted emission component might be expected if the out- 
flow arises in gas that is not a simple expanding shell. A conical 
outflow, for example, could produce a redshifted emission com- 
ponent and a blue-shifted absorption component, if viewed at 
orientations along the flow axis of the cone. The Fe line compo- 
nents from H-like Fe are consistent with an origin in an outflow 
of velocity of ~ lOOOkms -1 . The emission component of this 
H-like Fe pair is too narrow and weak to show up in the XMM 
pn spectra, but the pn spectra taken at approximately the same 
time and flux-state (the mid-state, Fig.[TJ are consistent with 
the presence of such a component. 

Having tight constraints on line energies and widths from 
HEG we tested the 2001 XMM pn and HETG data for consis- 
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Table 1. Spectral Lines detected in HETG data. Fixed parameters are denoted (f) 
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Fig. 3. X-ray spectra of NGC 3516 from XMM 2006 corre- 
sponding to the flux states denoted in Fi. 1, i.e. spectra are from 
the high/flare state (top), the mid-flux state (middle) and dip 
state (bottom). 



tency with the presence of H-like and He-like absorption lines 
at the same observed energies and widths as those measured 
during 2006. The pn spectra were inconsistent at the 1% signif- 
icance level with the presence of absorption lines at the same 
energies, widths and depths as observed in 2006. Specifically, 
the XMM data limit the presence of H-like Fe absorption to be 
an order of magnitude weaker during Nov 2001 and more than 
a factor 2 weaker during April 2001 than was observed during 
2006. Fig. [5] shows 2001 November pn spectral data overlaid 
with that from 2006, illustrating the difference in spectral sig- 
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Fig. 4. X-ray spectra of NGC 3516 from XMM 2001 April, 
2001 Nov, Suzaku data from 2005 October and 2006 Oct XMM 
data 



natures for the two epochs. We also analyzed the Suzaku data 
for NGC 3516. The data were reduce d following the analy- 
sis procedure of iMarkowitz et al.ll2007l and the XIS 0, 2 and 3 
spectra were summed for spectral fitting. The data were fit for 
the presence of the H-like Fe absorption line and fitting ruled 
out the presence of a line of the same flux or equivalent width as 
observed during 2006, at the 1 % significance level (assuming 
the line energy and width to be as detected by HETG). 

We then searched the low energy part of the 2006 HETG 
data, fitting MEG and HEG data together to search for nar- 
row features. Several lines are significantly detected (Fig.|7]i 
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Fig. 5. XMM pn data shown as a ratio to the local continuum 
fit in the Fe K regime during Oct 2006 (red) and Nov 2001 
(black). 

and were modelled with a simple gaussian form. Table 1 lists 
lines detected across ~ 1 - 8 keV in HETG data. Line ener- 
gies tabulated have been corrected for the systemic redshift 
(z=0. 008836). Photon fluxes are given with positive values 
indicating an emission line and negative an absorption line; 
equivalent widths are measured relative to the local observed 
continuum. The velocity implied for the gas is noted for each 
detected line, where negative velocity indicates outflow (no ve- 
locity is given for the 6.7 keV line as this is an uncertain blend 
of ionization states). The last column of Table 1 shows the x 1 
value when the line is added to the model. Quoted errors are 
90% confidence intervals. The Fe KJ3 line in Table 1 was scaled 
to the Ka component but was not significantly detected in its 
own right. The K/3 energy was fixed appropriate to an origin in 
neutral gas and constrained to have flux 1 1.3% that of the neu- 
tral Fe Kc line component. For groups of lines with closely- 
spaced energies we linked the line widths together in the fit so 
a turbulent velocity could be estimated based upon several lines 
together; for lines whose widths were linked the same value is 
shown for the <x entry for each but only the first occurrence in 
Table 1 carries the combined error. 

The fitted widths correspond to velocity dispersion 
o- Mg = 1245+^°kms _1 , a- Si = 1265+j^kms^ 1 and cr s = 
1292^25k ms ' : so ' together the lines indicate a velocity dis- 
persion <x ~ 1300 km s . The weighted mean outflow velocity 
from the absorption line measurements is ~ 1575 +254 kms~'. 

We examined the RGS data in the 1-2 keV band to seek 
confirmation of the features detected. RGS data confirm the 
presence of Mg xi and xn absorption lines (absence of each 
is rejected at p < 1% significance), yielding line energies 
E = 1358+°;™* E = 1.482 + 0.005 keV with correspond- 
ing bulk velocities 1239 and 2000 km s for the outflows. The 
data are consistent with the detection of the Si xm absorption 
line but RGS has so little area at this energy that no useful line 



constraints could be obtained. RGS has a slightly lower energy 
resolution than the HETG and could not detect the weak emis- 
sion lines in this regime. 

To put the new results into context, we now fit the broad- 
band pn spectrum utili zing models for the photoionized gas 
run from xstar v211n ( Kallman et all 2004 ). As this source is 
known to show a large degree of flux-related spectral variability 
we separated the pn data into high, mid and low-state spectra, 
as indicated in Fig. Q] For the initial fit we took the spectral data 
encompassed by the mid-flux state and fit across 1-10 ke V. Data 
in this band require inclusion of four layers of gas in addition 
to a small column of neutral gas fi xed at the Galactic line-of - 
sight value (N H = 4.08 x 10 2() cm' 2 , lDickev & Lockmanll990l) . 
Table 2 details the results of fitting the pn spectra covering the 
mid-level flux (note that the constraints on the ionization pa- 
rameter of zone 4 were unique in being derived from the HEG 
and not the pn fit). 

The two low-£ layers of gas show signatures in the RGS 
data below 1 keV and provide some of the broad spectral cur- 
vature. These two layers (Table 2: zones 1 & 2) show a similar 
total column density to that typically detected in the UV band 
for NGC 3516. No significant features are expected from zones 
1 or 2 above ~ 1 keV and as this gas is likely to exist relatively 
distant from the active nucleus we do not discuss those zones 
in detail in this paper. 

In the pn fit we found all layers of gas to be consistent with 
full covering of the source except for zone 3 (Table 2) whose 
covering fraction is 45 + 4% for gas with Nh ~ 2 x 10 23 cm -2 , 
log £ ~ 2. The most extreme conditions detected are in a 
layer producing the He-like and H-like Fe absorption lines 
where the line depths can be fit with a zone of column den- 
sity Nh = 2.62 x 10 23 cm~ 2 . Overall this 4-zone model gave 
X 1 = 2093/1746 d.o.f with most of the contribution to^ 2 re- 
maining in the soft-band near 2 keV. This may be indicative of 
a need to refine the model but as additional features indicated 
by the pn close to 2 keV do not show up in the contemporane- 
ous HETG data, and as this is also a region most prone to pn 
calibration uncertainties we do not attempt to further define the 
model based on that part of the spectrum. 

Having found a reasonable solution for the source it is 
interesting to revisit the individual lines detected in HETG. 
Examination of the xstar solution for zones 3 and 4 is partic- 
ularly interesting. Zone 4 is of such high ionization-state that 
little absorption arises from that gas below the Fe K regime. It 
is most likely that the lines detected in the 1-3 keV band arise 
in zone 3, the partially-covering absorber. We conclude from 
the combined HETG and pn results that the partially-covering 
zone 3 gas has an intermediate ionization-state and large col- 
umn composed of clouds with velocity widths cr ~ 1210+1210 
kms^ 1 and bulk outflow velocity ~ 1575 + 254 km s -1 . 

A model utilizing an ionized reflector (parameterized us- 
ing the reflion model of Ross & Fabian 2000) in place of 
zone 3, but allowing in the fit the other three layers of warm 
absorbers did not fit the mid-state spectrum at all well (y 2 = 
3441/1714 d.o.f.). While some reflection from the absorbing 
clouds must (at least) be present in this source it appears the 
absorption effects are dominating the spectral shape and any 
additional reflection was not discernible in these data. XMM 
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lacks the broad bandpass needed to separate out any reflection 
component from the absorbers and reflection-dominated mod- 
els are not considered further in this paper. 

3.1. Spectral Variability 

Havin g confirmed t he pre viously-established picture for NGC 
3516 (Tur ner et all 2005) as a source dominated by layers of 



Zone column density 

( 10 22 atoms cnT 2 ) 



log (^/ergcms ') 



complex absorbing gas, we then investigated the spectral vari- 
ability of the source in the context of such a model. 

First, we attempted to model the spectral steepening 
with flux by allowing each of the key absorber parameters 
(ionization-state, column density and covering fraction) for 
each layer to vary in turn. We found only one way in which a 
variation of a single fit parameter could adequately model spec- 
tra exhibited over the full flux range traversed by the source 
and that was a change in covering fraction of zone 3. As noted 
previously, the mid-state data represent a covering fraction of 
45 + 4% for zone 3. A good fit can be achieved for the high- 
state spectrum by allowing a simple drop in covering fraction 
to 37 + 5% and for the low state by requiring an increase to 
67 + 7% for zone 3. All other single-parameter explanations for 
the observed spectral variability could be ruled out at p < 0.1% 
significance (mostly failing to fit the low state). Specifically, the 
popular model where the relative strengths of a variable power- 
law and steady reflection component change cannot account 
for the deep dip spectrum (this is ruled out at p < 0.1% signifi- 
cance). In the context of the variable-coverin g model for spec 



tral v ariability then the Suzaku observation (Mar kowitz et al. 



20071) caught NGC 3516 hidden by a large cool cloud - the 
overall behavior of NGC 3516 appears to be a mixture of 
cloud passages with changes in ionization-state of the absorb- 
ing zones. What is perhaps most surprising is the rapidity of 
some of the changes observed. 

Using crude time-resolved spectroscopy we found that for 
both the flare and dip events, the change in covering fraction ac- 
counts for all but 15% of the source flux variations that are ob- 
served, leaving a question as to whether all observed variability 
can be caused by simple absorption changes. We therefore con- 
structed a softness ratio to investigate the spectral variability 
at finer timescales than allowed by time-resolved spectroscopy. 
Fig.[8]shows the softness ratio (0.5-4.0/4.0-1 0.0 keV countrate) 
versus the total 0.5-10.0 keV count rate compared to a model 
line representing the values allowed by our multi-zone model, 
for different values of covering fraction (assuming all other 
parameter values to be fixed, including intrinsic continuum 
normalization). The model line terminates at the point where 
the source is completely uncovered. The source behaviour is 
broadly consistent with the model line in the range range 10-15 
ct/s. Below 10 ct/s there is some residual soft flux representing 
some unmodeled soft emission. Above 15 ct/s there are some 
data showing the source to be softer than expected, these data 
may represent times when the ionization-state of one of the gas 
layers has increased (the ionization-state of absorbers laying 
outside of zone 3 would be expected to increase as the cover- 
ing fraction of zone 3 went down because of changes in the 
incident illuminating flux). Flux changes above 20 ct/s are con- 
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Table 2. X-ray Absorber Parameters and 90% confidence er- 
rors, fixed parameters are denoted (f) 



sistent with no spectral variation and may represent intrinsic 
continuum variability. 




6.5 7 
Observed Energy (keV) 



Fig. 6. Spectral data from the summed HEG first-order data 
from the five data segments taken during Oct 2006. The solid 
line shows the model fit to the data. 




.5 2 2.5 

Observed Energy (keV) 

Fig. 7. Summed first-order MEG data and summed first-order 
HEG data from the five segments accumulated during Oct 
2006. The two grating datasets are shown coadded, the solid 
line shows the model detailed in Table 1 . 

Finally, we note the presence of weak transient Fe emission 
features in the spectra, o f the type first discov ered in the 2001 
observation of this AGN (ITurner. et all 120021) . The addition of 
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a narrow (cr ~ 3eV) gaussian emission line at E — 6.20 + 
0.02 keV was required in the pn data (absence of a line being 
rejected at p < 1% significance), with line flux n = 9.4 + 1.4 x 
10~ 6 photons cm"~ 2 s -1 , EW=23+4 eV. Addition of a second line 
to the model yielded E = 5.98 + 0.03 keV, n = 6.0 + 2.0 x 
10~ 6 photons cm _2 s _1 , EW=15+5 eV. Some residuals are also 
evident in the HEG spectrum at those energies supporting the 
reality of the features (Fig.|6]l. Interestingly, the features show 
up most prominently in the high-state spectrum (Fig. [3), as they 
did during the 2002 observation. This tendency for Doppler- 
shifted Fe emission lines to show up in source high flux states 
supports a close link to the continuum source flux with little 
lag, ie an inner disk origin. 
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Fig. 8. The ratio 0.5-4.0 keV/4.0-10.0 keV pn cts on the y-axis 
versus the total pn count rate 0.5-10.0 keV on the x-axis. The 
solid model line shows the changes expected if the only vari- 
able is the covering fraction of zone 3 gas, the model line ter- 
minates at the point where none of the continuum source is 
covered by zone 3. 



4. Discussion 

Here we detect four distinct layers of ionized gas absorbing the 
nuclear continuum. The three zones detected during 2001 have 
obvious counterparts in 2006 data, and a new, highly ionized 
layer is detected for the first time in this source. In our new 
analysis we also find two low-ionization zones of gas, similar 
to those known to cause absorption in the UV regime, and as 
these likely arise quite far from the nucleus, are not the topic of 



detailed discussion here. We concentrate on the implications of 
the findings for the innermost zones, 3 and 4. 



4.1. The Partial-Covering Gas: Zone 3 



The 'h eavy' partial-covering zone detected by iTurner et al 



(2005) is very similar in all respects to the layer denoted "zone 
3" here: this is a zone characterized by intermediate ionization- 
state (log £ ~ 2), large column density (Nh ~ 2 x 10 23 crrT 2 ) 



and modest covering fraction, varying around 
to that observed in 2001. 



45%, similar 



The source light curve shows a deep dip in flux having an 
egress of duration ~ 30 ks. Time resolved spectroscopy indi- 
cates that the simplest explanation of the spectral variability 
observed during the dip is a simple change in covering fraction 
of the zone 3 absorber. In the context of this model the Suzaku 
spectra may be interpreted as a long duration eclipse by a cool 
gas zone, leading to an associati on of the Suzaku 'primary 
absorber' dMarkowitz et all 120071) with zone 3. Examination 
of the softness-flux plot (Fig. [8} indicates that while covering 
changes in zone 3 may cause deep observed dips in X-ray flux, 
other important phenomena are probably at work to produce 
the total light curve and spectral variations observed in the 
source. Possibilities include ionization changes of the absorber 
and some intrinsic variation of the source continuum flux. 

The evidence for rapid changes in covering fraction of 
an X-ray absorber places this BLSyl galaxy in the com- 
pany of sources such as the narrow-line Seyfe r t 1 (N LSyl) 



galaxi es MCG-6-3 0-15 ( McKernan & Yaq oobi Il998l) . Mrk 
766 dTurner et all [2007b, and the BLSyl 1H04 19-577 
dPounds et all 20041) : all of which have been discussed by some 
authors in terms of variable absorption. Partial-covering mod- 
els have been raised at various epochs as a possibl e explana- 
tion f or Seyfert spectral and timing properties (e.g. lBoller et al 
19971) . the complete ambiguity between absorption and reflec- 
tion dominated models in low-resolution X-ray spectra avail- 
able to date left it unclear which process was most dominant. 
However, recent results have shown that the role of absorption 
has not been adequately accounted for and consideration of 
rapid absorption variability is obviously key to any meaning- 
ful analysis of the timing and spectral variability properties of 
active galactic nuclei. 

The presence of large amounts of optically-thick gas close 
to the active nucleus means that both absorption and reflection 
from that material must contribute to the observed spectrum 
and in isolating key events, careful decomposition of the data 
train may lead to significant new insight into the nature of the 
reprocessed emission. 

As noted previously, the detailed shape of the dip event in 
NGC 3516 is complex (Fig.©. Some fraction of the source is 
completely covered for 10 ks and then the egress begins with a 
20 ks shallow recovery followed by a rapid rise over 6 ks. The 
shape of the dip is striki ngly similar to that observe d for a deep 
dip in MCG-6-30-15 dMcKernan & Yaqoobi [19981) . indicating 
that the same physical process is likely responsible for deep 
flux dips across the Seyfert population. 
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If one considers a spherical cloud of gas of constant den- 
sity occulting a point continuum source, although the cloud 
shape produces some gradation in opacity as the blob traverses 
the sight-line, this crossing would produce a smoothly varying 
light curve, and symmetric in ingress/egress times, rather than 
the asymmetric curve that is observed. 



AR 9.2 x 10 12 v, 4 cm and AR/R 2.6 x 10 



-3,, -2 



indicating 



As suggested for MCG-6-30-15 by iMcKernan & Yaqoob 
dl998h the dip profile shape can be explained by a contin- 
uum source having non-uniform emissivity. If the emitter is 
composed of a central bright source embedded in a fainter 
extended emission region, one could get a dip profile sim- 
ilar to that observed as a cloud crosses the line-of-sight. 
Mc Kernan & Yaqoobl ( 1998b sketch a schematic of the gen- 



eral picture (their Fig. 2), applicable also to NGC 3516 (al- 
though we require no gap between the central and extended 
components of emission). The dip profile observed in NGC 
3516 is consistent with the IMcKernan & Yaqoobl dl998l) sce- 
nario. Another possibility is that inhomogeneities exist in the 
absorber, and that this complex absorber crosses a simple point- 
like continuum source. The asymmetry of the ingress/egress 
timescales gives us some additional information: these can ar- 
guably be more naturally produced with an inhomogeneous ab- 
sorber crossing the sight line than with a simple cloud crossing 
a spatially-complex emitter. In the latter case we would need 
to invoke additional complexity such as that the central source 
is asymmetric and that the occulting material does not simply 
cross the source but maybe changes trajectory during its pas- 
sage across the line-of-sight egressing across a different spatial 
dimension of the continuum emitter than it ingressed across. 
Complex situations such as this might arise if a cloud of gas is 
lifted up off the disk as it orbits, ie if the absorber is part of a 
disk wind. 

HETG data trace individual lines arising in zone 3, showing 
a mean outflow velocity ~ 1600 km s _1 , similar to the outflow 
velocity ~ 1 100 km s _1 found for what appears to be the same 
zone observed during 2001. If this zone appears to be occulting 
the central emission/reflection regions then the gas must have 
radial and transverse components of motion. 

As the transverse component of velocity is unknown, we 
express our derived quantities in units likely applicable to such 
sources. Using a mass estimate for t he central black hole, M= 
2.95 x 10 7 M o dNikolajuk et alll2006h we obtain an estimate of 
the radial location of zone 3, R - 3.8 x 10 15 v 4 2 cm where V4 is 
the orbital velocity in units of 10 4 km s _1 . For a velocity of 10 4 
km s _1 R ~ 1.5 l ight days, locating th e material on the inner 



edge of the BLR dWanders et all 1 1993) 



The fitted ionization parameter, £ may be related to the ra- 
dial location of the gas zone, R, the gas density n and the (in- 
trinsic) ionizing luminosity using £ = Li on /R 2 n and thus can be 
used to estimate n (these, and other simple relationships used 
here are derived in detail by ICrenshaw. Kraemer & George! 
2003 and lBlustin et al.ll2004l) . The ionizing luminosity utilized 
by xstar is defined over the l-1000Ryd band. The observed 
luminosity is L Q b s - 1.2 x 10 43 ergs~'; correcting for the 
fitted absorption layers gives an implied intrinsic luminosity 
Li on =i 4.8 x 10 43 erg s _1 . Utilization of the implied intrinsic lu- 
minosity gives n ^ 2.1 x 10 10 v^crrT 3 . Combining this with the 
relation Nh — nAR where AR is the shell thickness, we obtain 



zone 3 to be a thin shell of gas. 

For zone 3, the estimate of radial location and comparison 
of assumed transverse and measured outflow velocity indicates 
that the gas could be gravitationally bound. While a transverse 
velocity of order ~ 10 4 km s _1 is necessary to obtain a reason- 
able geometry, the value is subject to large uncertainty. Further, 
as we only measure the outflow component along our line-of- 
sight the true outflow velocity may be larger than observed and 
zone 3 is thus consistent with having comparable transverse 
and outflow velocities. In addition to this geometrical consid- 
eration we note that steady-state wind model solutions exist 
where there is continuous injection or conversion of momen- 
tum into the wind, and it is not strictly necessary to have a 
measured outflow velocity that is greater than the local escape 
velocity for an outflow ultimately to be unbound. Hence we 
cannot at this stage make a definite statement about whether the 
observed outflow forms part of a large-scale wind or whether it 
indicates a more complex gravitationally-bound flow. 

The mass in the absorber can be estimated from M 
AnR 2 NHm P C g where m p is the mass of a proton and C g is 
the global covering fraction. The mass in zone 3 is estimated 
as 0.03 v^ 4 C,M n . The mass flow rate is then M M/t c 
dCrenshaw. Kraemer & Geor ge. 2003) where t c - AR/v is the 



radial travel time and here v is the outflow velocity rather than 
the orbital velocity. As noted above, for zone 3 we have mea- 
sured an outflow velocity of 1600 km s _1 that is is a lower limit 
on the actual outflow velocity, as it is only the component of 
velocity that is in our line-of-sight. We estimate the mass flow 
rate as 16vi6oo C g M year -1 , this is quantity that does not 
depend on the unknown orbital velocity. 

Assuming the dip event was caused by a simple cloud, its 
diameter d c i ol ,d — v or t,t, where t is taken as the time for the 
egress, the strongest constraint within the dip event, so d c i m ,d 
3xl0 13 V4 cm. For a velocity of 10 4 km s _I the lateral dimension 
of the cloud (d c i ou( i) is then approximately equal to the cloud 
depth AR and a reasonable cloud geometry is thus consistent 
with the likely orbital velocity. The maximum coverage by the 
absorber is ~ 70% suggesting the continuum to be ~ 1 .2 times 
the size of the cloud on the sky. The light-crossing time for the 
continuum source is thus ~ 10 3 V4 s, consistent with observed 
variability. 

4.2. The Highly Ionized Gas: Zone 4 

The highly-ionized zone 4 gas detected here has log £ ~ 4.3 
and column density Nh ~ 2.6 x 10 23 crrT 2 , adding to the evi- 
dence for large columns of highly-ionized gas being important 
in Seyfert nuclei. This gas zone is the origin of the H-like and 
He-like absorption lines from the K shell of Fe but has little 
signature below ~ 6 keV. The 2001 XMM data are inconsistent 
with the presence of He-like and H-like absorption lines of the 
same flux. 

The very high ionization and column density, combined 
with apparent variability on a timescale of years and the en- 
ergies at which these lines are detected disfavors an origin of 
this gas in the warm-hot intergalactic medium. That we have 
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detected similarly extreme z ones of gas in the longest ob- 
servations of AGN (Mrk 766, iTurner et al.l 120071: NG C 3783, 
Reeves et al.ll2004t MCG-6-30-15. lYoung et all2005l) suggests 
their existence could be common in active nuclei, these fea- 
tures are below the limit of detectability for typical short AGN 
exposures conducted to date. 

An intriguing aspect of this extreme zone is the appar- 
ent observation of both the emission and absorption compo- 
nents from H-like Fe in zone 4. The pair of features appears 
reminiscent of classic 'P Cygni' profiles that arise in outflow- 
ing shells of gas such as stellar winds. In this case it is not 
clear that the components arise from such a simple spherical 
shell of gas. Most likely, the gas has some more complex ge- 
ometry (such as a bi-cone). Previous claims of 'associated' 
emission/absorption features have been made for the lower- 
ionization gas zo nes traced by soft-band X-ray features (e.g. 
Kaspi et al.l2002r) . From the velocity separation of the emission 



and absorption components of H-like Fe the implied outflow 
velocity of zone 4 is ~ 1000 kms 1 while line widths imply a 
velocity dispersion <x ~ 2124+?J!? kms -1 (based on the isolated 
H-like absorption line). 

Using the equation for ionization parameter as before we 
estimate a density for zone 4, n 1 .6 x 10 8 v\ 4 cirT 3 and AR 
1.6 x 10 15 v 4 ~ 4 cm, AR/R =s 0.4 v^ 2 . The mass in zone 4 is M ~ 
0.04v4 ~ 4 C g M . The measured outflow velocity for zone 4 is ~ 
1000 km/s and the mass loss rate 0.06viooo C g M year~'. 

Turning to the emission signature of zone 4, the normaliza- 
tion of the H-like Fe emission line can be compared to the emis- 
sion predicted by xstar for a full shell of gas covering the en- 
tire ionizing source. This comparison shows the observed line 
flux to be an order of magnitude larger than predicted. Taken 
at face value this leads us to an estimate of the global cover- 
ing factor (covering factor seen over An steradians as viewed 
by the continuum source) to be C g ~ H^g 7 for zone 4, which 
is unreasonably large, especially given the mass and loss rates 
estimated above. 



Despite the sizable uncertainty on this covering factor 
it is interesting to consider the ways in which the pre- 
dicted emission line strength may be brought into line with 
that observed. Since the Fe xxvi emission is a resonance 
line, its strength could be increased by photo-excitation, 
as suggested for the resonance lines of H and He-l i ke O, 
N, and Ne in several other Seyferts (e.g. ISako et. al 2000; 



Armentrout, Kraemer, & Turner 120071) . However, there will be 
little enhancement by photo-excitation once the line becomes 
optically thick to the continuum radiation. Although there can 
be greater contribution from photo-excitation if the gas is tur- 
bulent, the ionic column densities are too large in this case 
(~ 10 19 cirT 2 ). For example, if the absorbing gas has a micro- 
turbulence <x ~ 1500 km s _1 , then the predicted line flux would 
be twice as large as that calculated in the xstar table used here. 
However, this is not enough of an increase to bring the pre- 
dicted and observed fluxes into consistency. Even allowing for 
the absorption-correction to the intrinsic luminosity incident 
upon zone 4 doesn't produce any significant increase in the 
predicted line flux. The recombination line depends critically 
on the fractions of Fe xxvi and fully ionized Fe, and there is a 



limit to how much recombination line can be produced in the 
gas. 

Nevertheless, the strength of the Fe xxvi line is very sensi- 
tive to temperature, particularly since it is primarily formed by 
recombination. Our models predict that when that is the dom- 
inant Fe ion, the gas must be in the Compton-cooled regime 
dKrolik. McKee. & Tarter 1119811) . In this case, the gas temper- 
ature depends mostly on the spectral energy distribution of 
the incident radiation. A weaker X-ray flux or a stronger low 
energy flux could result in a significantly lower temperature, 
hence a stronger Fe xxvi line. For example, the presence of a 
reserve of low energy photons, such as may occur close to the 
inner edge of the accretion disk, can drive the Compton temper- 
ature down. In conclusion we note that in NGC 3516 Compton 
cooling appeai-s to be important. 

5. Conclusions 

Emission and absorption features are detected across the X- 
ray spectrum of NGC 3516 confirming the presence of at least 
four distinguishable layers of absorption that shape the spec- 
trum. Emission and absorption lines are detected from H-like 
Fe revealing the presence of one absorbing layer with extreme 



properties, log £ ~ 4.3, Nh ~ 2.6 x 10 cm outflowing at 
~ 1000 km s" 1 . 

A lower-ionization zone of gas with log £ ~ 2.2, Nh ~ 
2xl0 23 cirT 2 covers just ~ 45% of the source. Variations in cov- 
ering fraction of this intermediate layer account for much of the 
spectral and observed flux variability in the source. In particu- 
lar, a deep dip in the light curve can be explained as an eclipse 
by a cloud of this gas and the shape of the light curve during 
occultation reveals that either the cloud or the emission region 
must be inhomogeneous. The remarkable similarities between 
the deep minima of NGC 3516 and MCG-6-30-15 suggest that 
these events might be common in Seyfert galaxies. It appears 
that rapid variations in absorber parameters play a much greater 
role in shaping the observed X-ray variability in AGN than has 
been appreciated to date. These occultation events may offer 
the best way forward for probing the details of the innermost 
regions of AGN. 

The detection of a significant outflow velocity for the most 
highly-ionized gas zone in NGC 3516 suggests that at least this 
component of X-ray absorption originates in a wind. The small 
implied radial location for the highly-ionized gas rules out the 
putative torus as an origin fo r the wind, leavin g a disk wind 
as the most likely explanation. Ruiz et al.l ( 2005 ) suggest, from 
optical kinematic studies, that NGC 3516 appears to be viewed 
at a fairly high inclination, ie close to the plane of the absorb- 
ing material and accr etion disk system, similar to NGC 4151 
(ICrenshaw et al. , 2000). Sources viewed at shallow angles to 
the plane of the disk would be the best candidates in which 
to view disk winds, and this may have contributed to the easy 
detection of outflows in highly-inclined sources such as these. 
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